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The  Effects  of  a  Low-Altitude  Nuclear 
Burst  on  Millimeter  Wave  Propagation 


I.  INTRODUCTION 


In  this  report  the  limitations  imposed  on  millimeter  waves  by  a  low  altitude 
nuclear  burst  are  examined.  Associated  with  the  burst  is  a  fireball,  an  extremely 
hot  and  highly  ionized  spherical  mass  of  air  and  gaseous  weapon  residues.  The 
fireball  grows  rapidly  and  because  of  its  intense  heat  some  of  the  soil  and  other 
material  in  the  area  are  vaporized  and  taken  into  the  fireball;  strong  afterwinds 
cause  large  amounts  of  dirt  and  debris  to  be  sucked  up  as  the  fireball  rises. 

In  order  to  estimate  the  effects  of  dust  on  the  performance  of  a  millimeter 
wave  system  it  is  necessary  to  first  develop  a  dust  model.  Since  there  is  a  great 
deal  of  uncertainty  as  to  how  representative  the  model  is  of  the  true  nuclear  environ 
ment,  a  sensitivity  analysis  of  the  attenuation  dependence  on  pertinent  dust 
parameters  is  conducted.  It  is  assumed  that  the  dust  particle  sizes  follow  a  power 
law  distribution.  Then  the  attenuation  is  computed  as  a  function  of  minimum  par¬ 
ticle  radius,  maximum  particle  radius,  particle  size  distribution,  index  of  refrac¬ 
tion  and  density.  A  Mie  formulation  is  used  to  calculate  the  absorption  and  scatter¬ 
ing  losses  produced  by  dust. 

The  total  attenuation  produced  by  a  1  Megaton  surface  burst  is  then  computed 
using  the  WESCOM  code.  The  attenuation  consists  of  contributions  due  to  fireball 
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ionization,  dust  and  atmospheric  oxygen  and  water  vapor.  Attenuations  are  obtained 
as  a  function  of  time  after  burst,  distance  from  burst,  elevation  angle  and  fre¬ 
quency. 


Z  PARTICLE  SIZE  DISTRIBUTION 

The  sizes  of  dust  particles  are  often  represented  by  a  power  law  probability 
distribution  of  the  form 

P(r)  =  kr”p  (1) 

where  r  is  the  particle  radius,  p  is  the  exponent  and  k  is  a  constant,  selected  such 
that  P(r)  is  a  proper  probability  distribution.  For  some  applications  a  log-normal 
distribution  of  particle  sizes  is  used  since  it  provides  a  better  model  of  the  very 
small  particles.  For  this  study,  however,  we  will  show  that  the  actual  distribution 
of  the  very  small  dust  particles  does  not  significantly  affect  the  attenuation  so  the 
power  law  distribution  is  suitable. 


dr  =  1  . 


Solving  for  k  we  have 


r-lp-l>  _  -ip-1) 
min  max 


where  rm^n  and  rmax  are  the  minimum  and  maximum  particle  radii  respectively. 
The  total  number  of  particles  of  radius  r  is  then 


N(r)  =  NTP(r) 


where  N„  is  the  total  number  of  dust  particles. 


1  •  *  v  ' >  »  w  * 


and 


Pb  =  bulk  density  of  particulate  material, 
V  =  mean  volume  of  a  particulate. 


U 


Then 


N(r)  = 


pb  (4/3  irr11) 


_ __  r'P 


(p-1)  _  -(p-1) 


min 


max 


(6) 


f. 

Cr' 


It  can  be  shown  that 


r-(p-4)  .  r-(p-4) 

3  p-1  min  max 

p7*  '  “-(p-ir;  r-(p-i)  • 
min  max 


Finally, 

N(r)  = 


Mt  (p-4)  r 


-P 


4/371  p. 

"pr'4>  7 

r-(P“4)1 

b 

[  mm 

max  j 

(8) 


3.  ATTENUATION  OF  DUST  PARTICLES 

Millimeter  waves  incident  on  atmospheric  particulates  undergo  absorption  and 
scattering,  the  degree  of  each  being  dependent  on  the  size,  shape  and  complex 
dielectric  constant  of  the  particle  and  the  wavelength  and  polarization  of  the  wave. 
An  expression  for  calculating  the  absorption  and  scattering  from  a  dielectric  sphere 
was  first  derived  by  Mie.  1  It  has  the  form 


r: 

L' 


Q.  =  -4r  Re  E  (2n  +  1)  (as  +  bS) 

t  2tt  n=i  n  n 


(9) 
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1.  Mie,  G.  (1908)  Contribution  to  the  optics  of  suspended  media  specifically  colloial 
metal  suspensions  (in  German),  Ann.  Physik,  ^25:377  -445. 
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where  Q  ,  the  extinction  cross  section,  represents  losses  due  to  both  absorption  Q  . 

and  scat+ering  Q  ,  and  a  and  b  are  very  complicated  functions  of  the  spherical 

Bessel  tc:ms  uiat  correspond  to  the  magnetic  and  electric  modes  of  the  particle 

2 

respectively.  Q,  has  the  dimension  of  area  and  is  usually  expressed  as  cm  . 

1  2 

Physically,  if  a  wave  having  a  flux  density  of  S  W/cm  is  incident  on  the  particle, 
then  S  X  is  the  power  absorbed  and  scattered. 

When  the  size  of  a  dust  particle  is  very  small  with  respect  to  wavelength,  then 
the  Rayleigh  approximation  is  valid.  For  this  case 


(10) 


Qs  = 


128  jt5  r6 
3~X^ 


n2-l 

nK  2 


(11) 


It  is  seen  that  the  absorption  is  inversely  proportional  to  the  wavelength  while  the 
scattering  loss  is  inversely  proportional  to  the  fourth  power  of  the  v/avelength;  thus 
when  the  wavelength  is  large  compared  to  the  particle  size,  the  absorption  domin¬ 
ates  and  scattering  losses  are  often  assumed  negligible.  Furthermore,  since  the 
absorption  is  proportional  to  the  particle  volume,  the  total  attenuation  is  propor¬ 
tional  to  the  total  volume  of  dust. 

As  the  dust  particles  become  larger,  then  the  Rayleigh  approximation  is  no 

longer  valid  and  the  Mie  formulation  must  be  used.  Both  absorption  and  scattering 

cross  sections  continue  to  increase  with  particle  size.  Finally,  after  reaching  a 

peak  the  total  cross  section  begins  to  level  off  and  would  eventually  approach  a 

value  of  twice  the  geometric  cross  section  of  the  particle  when  it  is  very  large  with 

2 

respect  to  wavelength.  Thus  we  note  that  as  the  particle  becomes  larger  the 
extinction  cross  section  which  was  initially  proportional  to  the  particle  volume 
becomes  proportional  to  the  particle  cross  sectional  area. 

The  attenuation  coefficient  is  equal  to 

oo 

a  =  J  N(r)  Q  dr  .  (12) 

o 

_3 

This  expression  is  in  nepers /cm;  if  r  and  Q,  have  units  of  cm  and  N(r)  is  in  cm  . 

*  5 

In  order  to  convert  to  dB/km  a  multiplicative  factor  of  4.  343  X  10  must  be 
introduced. 


2.  VanDeHulst,  H.  C.  (1957)  Light  Scattering  by  Small  Particles,  Wiley,  New  York 
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Q 

Attenuation  and  albedo,  the  ratio  of  scattered  loss  to  total  loss  ,  for  a 

Wt 

particle  size  power  law  distribution  are  plotted  in  Figure  1  for  sand  and  clay 
for  frequencies  from  10  to  95  GHz,  It  is  assumed  that  the  minimum  and  maximum 
particle  radii  are  0,  005  and  5  mm  respectively,  that  the  pov/o*-  lev  exponent  is 

3 

-3.  5,  and  that  the  densities  of  the  dust  material  and  average  mass  are  2.  6  gm/cm 

3 

and  100  gm/m  respectively.  It  is  seen  that  the  attenuation  increases  very  rapidly 
with  increasing  frequency.  It  should  be  pointed  out  that  attenuation  is  plotted  on  a 
logarithmic  scale  so  the  apparent  linear  increase  in  attenuation  at  the  higher  fre¬ 
quencies  actually  corresponds  to  an  exponential  increase,  in  reality.  Since  clay 
has  a  higher  index  refraction  than  sand,  the  attenuations  are  significantly  higher. 


Figure  1.  Attenuation  and  Albedo  of  Sand  and  Clay  Particles  at 
Millimeter  Wavelengths 


The  ratio  of  the  scattering  losses  to  the  total  losses,  often  referred  to  as  the 
albedo,  is  also  plotted  in  Figure  1.  It  is  interesting  to  note  that  the  albedo  for 
sand  is  approximately  0.  7,  even  for  a  wavelength  as  long  as  30  mm  when  the  aver¬ 
age  particle  radius  is  only  0.  02675  mm,  less  than  0.  001  wavelength.  This  indi¬ 
cates  that  even  though  the  average  dust  particles  are  very  small  with  respect  to 
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wavelength  the  scattering  losses  are  certainly  not  negligible  as  compared  to  the 
abso  "ption  because  of  the  limited  number  of  large  particles  that  are  present.  This 
issue  is  discussed  later  when  the  results  of  the  albedo  dependence  on  maximum 
particle  radius  are  presented. 

3.1  Dependence  of  Dust  Attenuation  on  Minimum  Particle  Radius 

In  Figures  2a  and  2b,  the  attenuations  are  plotted  for  sand  and  clay  respectively 
as  a  function  of  minimum  particle  radius  for  a  set  of  frequencies  from  10  to 
95  GHz,  It  is  seen  that  the  attenuation  rises  only  very  slightly  as  the  minimum 
particle  radius  is  increased.  It  is  interesting  to  note  that  it  passes  through  a  very 
broad  peak  at  95  GHz  and  would  probably  behave  similarly  at  the  lower  frequencies 
for  larger  minimum  particle  radii.  On  the  basis  of  these  results  it  can  be  con¬ 
cluded  that  the  attenuation  is  not  very  sensitive  to  the  minimum  particle  radius. 


Figure  2a.  Attenuation  at  Millimeter  Wavelengths 
as  a  Function  of  Minimum  Particle  Radius,  Sand 


6 


MINIMUM  PARTICLE  RADIUS  (mm) 

Figure  2b.  Attenuation  at  Millimeter  Wavelengths 
as  a  Function  of  Minimum  Particle  Radius,  Clay 


3.2  Dependence  of  Dust  Attenuation  and  Albedo  on  Maximum  Particle  Radius 

In  Figures  3a  and  3b,  the  attenuations  and  albedos  are  plotted  for  sand  and 
clay  respectively  as  a  function  of  maximum  particle  radius  for  the  same  set  of 
frequencies  from  10  to  95  GHz.  It  is  seen  that  the  attenuation  and  albedo  both  rise 
very  abruptly  as  the  maximum  particle  radius  is  increased.  The  attenuation 
reaches  a  peak  and  the  albedo  begins  to  level  off  when  the  maximum  particle 
diameter  is  approximately  equal  to  the  wavelength.  As  mentioned  previously  it  is 
'  .en  that  the  scattering  losses  are  comparable  and  in  some  instances  even  larger 
than  the  absorption  losses.  This  is  due  to  the  fact  that  the  imaginary  component 
of  the  index  of  refraction  is  very  small,  thus  resulting  in  a  small  absorption. 
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MAXIMUM  PARTICLE  RADIUS  (mm) 


Figure  3a.  Attenuation  and  Albedo  at  Millimeter  Wavelengths  as 
a  Function  of  Maximum  Particle  Radius,  Sand 


MAXIMUM  PARTICLE  RADIUS  (mm) 


Figure  3b.  Attenuation  and  Albedo  at  Millimeter  Wavelengths  as 
a  Function  of  Maximum  Particle  Radius,  Clay 
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ALBEDO 


3  J  Dtpndtact  of  Dut  Atteaaation  oa  the  Power  Law  Exponent 

In  Figures  4a  and  4b,  attenuation  is  plotted  as  a  function  of  power  law  exponent 
for  sand  and  clay  respectively  for  the  set  of  frequencies  from  10  to  95  GHz.  When 
the  exponent  p,  is  small  the  probability  of  having  large  particles  in  the  distribution 
is  relatively  high.  As  p  becomes  larger  fewer  large  particles  are  present.  Thus 
the  attenuation  passes  through  a  very  broad  peak  for  low  values  of  p  and  then  drops 
off  very  rapidly  when  p  becomes  larger.  Thus  the  power  law  exponent  can  sig¬ 
nificantly  affect  the  attenuation. 


Figure  4a.  Attenuation  at  Millimeter  Wavelengths 
as  a  Function  of  Power  Law  Exponent,  Sand 
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Figure  4b.  Attenuation  at  Millimeter  Wavelengths 
as  a  Function  of  Power  Law  Exponent,  Clay 


3.4  Dependence  of  Dust  Attenuation  and  Albedo  on  the  Index  of  Refraction 

In  Figure  5,  the  attenuation  and  albedo  are  plotted  as  a  function  of  index  of 
refraction  n=n^  -  jn2,  for  the  set  of  frequencies  from  10  to  95  GHz.  Both  the  real 
and  imaginary  components  of  the  index  of  refraction  are  increased  linearly  so  that 
the  particulates  have  an  index  of  refraction  that  corresponds  to  dry  sand  for  the 
lowest  values  and  water  for  the  highest  values.  It  is  seen  that  the  attenuation 
increases  and  the  albedo  decreases  sharply  as  the  particulates  change  from  dry 
sand  to  clay  and  then  both  essentially  level  off  as  the  real  and  imaginary  compo¬ 
nents  are  increased  further.  The  attenuation  and  albedo  are  examined  more  closely 
in  Figure  6.  In  Figure  6a,  the  imaginary  component  is  fixed  at  n2=  0.  2  and  the 
real  component  is  varied  from  1.  6  to  3.  4;  in  Figure  6b,  the  real  component  is 
fixed  at  n^  2.  5,  and  the  imaginary  component  is  varied  from  zero  to  0.  6.  For 
both  cases  the  attenuation  increases  gradually  as  either  of  the  components  of  index 
of  refraction  is  increased.  Based  on  the  curves  for  the  albedo  it  is  seen  that 
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increasing  the  imaginary  component,  increases  the  absorptive  losses  while  in 
creasing  the  real  component  results  in  an  increase  in  the  scattering  losses.  ' 
it  can  be  concluded  that  the  attenuation  is  very  dependent  on  both  the  real  and 
imaginary  components  of  the  index  of  refraction  over  the  ranges  from  about 
2  <  n1  <  4  and  0  <  n^  <  ■  8. 


Figure  5.  Attenuation  and  Albedo  at  Millimeter 
Wavelengths  as  a  Function  of  Index  of  Refraction 
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Figure  6a.*,  Attenuation  and  Albedo  at  Millimeter  Wave¬ 
lengths  as  a  Function  of  Real  and  Imaginary  Components 
of  Index  of  Refraction.  Real  Component 
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Figure  6b.  Attenuation  and 
Albedo  at  Millimeter  Wave¬ 
lengths  as  a  Function  of 
Real  and  Imaginary  Compo¬ 
nents  of  Index  of  Refraction, 
Imaginary  Component 


3.5  Dependence  of  Dnit  Attenuation  on  Dcnilty  of  Particlet 

In  Eq.  (5)  it  is  seen  that  the  total  number  of  dust  particles  is  directly  propor¬ 
tional  to  the  total  mass  of  the  particles  M^,,  divided  by  the  bulk  density  of  the 
individual  particles.  In  this  formulation  a  single  scattering  model  is  used  so  the 
attenuation  is  directly  proportional  to  the  number  of  particles  as  seen  in  Eq.  (12) 
and  therefore  directly  proportional  to  M^/p^.  Thus  any  uncertainty  in  either  the 
total  mass  or  particle  density  will  produce  a  corresponding  error  in  the  attenuation. 


4.  ATTENUATION  PRODUCED  BY  A  NUCLEAR  BURST 


As  mentioned  previously  associated  with  the  nuclear  burst  is  a  fireball.  The 
initial  electron  density  within  the  fireball  is  extremely  high  and  for  a  period  of 
about  20  sec  blacks  out  the  propagated  wave.  After  that  time,  beta  radiation  from 
the  radioactive  debris  within  the  fireball  may  sustain  sufficiently  high  ionization 
levels  to  absorb  millimeter  wave  signals  for  several  minutes.  Thus  the  two  prin¬ 
cipal  sources  that  attenuate  the  propagated  wave  are  fireball  ionization  and  dust. 

The  nuclear  explosion  dust  model  is  often  divided  into  five  distinct  regions 
which  are  modeled  somewhat  crudely  in  Figure  7.  The  ejecta  region  consists  of 
relatively  large  particles  of  dust  and  debris  thrown  out  of  the  crater;  settling 
occurs  quickly  and  its  importance  diminishes  in  about  a  minute  or  so.  The  blast 
wave  produces  a  low  level  dust  region  often  referred  to  as  the  pedestal  or  sweep- 
up  layer.  This  region  has  high  dust  densities  for  about  a  minute  also;  lower 
densities  exist  for  many  minutes  later.  The  cylindrical  stem  region  forms  after 
the  establishment  of  the  fireball  vortex  and  lasts  for  several  minutes.  The  mush¬ 
room  shaped  cloud  or  main  cloud  is  the  major  dust  region  and  exists  for  many 
minutes.  Finally,  the  main  cloud  transcends  into  the  fallout  region  which  exists 
for  a  long  period  of  time. 

In  this  study  the  Weapons  Effects  on  Satellite  Communications  code  (WESCOM) 
is  used.  It  utilizes  environment,  propagation  and  signal  processing  models 
developed  by  the  Defense  Nuclear  Agency  (DNA)  and  attempts  to  provide  a  best 
estimate  of  the  quantities  being  modeled.  Included  along  with  ionization  and  dust 
losses  are  the  oxygen  and  water  vapor  losses  of  the  ambient  troposphere. 

The  WESCOM  dust  model  is  a  somewhat  simplified  model  in  that  it  assumes 

that  all  the  dust  is  initially  contained  within  the  fireball  and  then  forms  a  main  cloud 

3 

which  eventually  becomes  the  fallout  region.  The  stem  and  pedestal  are  not  in¬ 
cluded  in  the  WESCOM  dust  model.  Dust  particulates  within  the  fireball  are  divided 
into  the  eight  parties  size  ranges  shown  in  Table  1. 


3.  Thompson,  .  .  (1980)  Dust  Clouds  -  Models  and  Propagation  Effects,  Pro¬ 
ceedings  of  Submillimeter  Atmospheric  Propagation  Applicable  to  fiadar  and 
Missile  Systems,  ftedstone  Arsenal,  Alabama^  TK-8U-3,  pp.  1 14-117. 
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Figure  7.  Dust  Regions 


Table  1.  Particle  Size  Distribution 


Dust 

Group 

Minimum 

Particle  Diameter 

Fraction 
of  Particles 

1 

0.  01  mm 

1-1.  372  X  10~6 

2 

0.  9 

1. 356  X  10'6 

3 

4.  0 

1.463  X  10'8 

4 

10.  0 

8.  750  X  10'10 

5 

20.  0 

9. 587  X  10"11 

6 

32.5 

2.  222  X  10' 11 

7 

52.  5 

4.540  X  10"12 

8 

75 

1.  370  X  10"12 

Each  group  of  particle  sizes  are  initially  assumed  to  be  uniformly  distributed 
throughout  a  disc-like  cylinder.  These  disc3  are  stacked  within  the  fireball  as 
shown  in  Figure  8a  with  the  largest  particles  contained  within  the  lowest  disc  and 
with  decreasingly  smaller  particles  in  the  upper  discs.  Each  particle  size  group 
rises  to  a  maximum  altitude  at  which  point  the  particles  start  to  fall  with  their 


respective  terminal  velocities.  A  simplified  model  of  the  size  and  height  of  the 
rising  fireball  as  a  function  of  time  is  shown  in  Figure  8b.  The  terminal  velocities 
are  plotted  in  Figure  9.  Thus  we  have  a  model  for  which  the  largest  particles 
remain  aloft  for  a  relatively  short  period  of  time  while  the  smaller  particles  may 
remain  suspended  for  many  minutes. 


Figure  8a.  Dust  Model, 
Disc- like  Regions  Within 
Fireball 


Figure  8b.  Dust  Model,  Rising  Fireball  as  a  Function  of  Time 
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Figure  9.  Terminal  Velocity  of  Dust  as  a  Function 
of  Particle  Radius 


A  set  of  attenuations  were  computed,  using  the  WESCOM  code,  as  a  function 
of  frequency,  elevation  angle,  distance  from  the  burst  and  time  after  burst  for  a 
megaton  surface  burst.  In  Figure  10  the  attenuation  is  plotted  as  a  function  of 
time  for  a  slant  path  30°  above  the  horizon,  distance  from  burst  of  10  km  and  fre¬ 
quencies  of  10,  20,  45  and  95  GHz.  It  is  seen  that  the  attenuation  has  essentially 
two  peaks;  the  first  peak  occurs  about  35  sec  after  the  burst  and  is  caused  by  the 
dust  rising  through  the  antenna  beam,  the  broader  and  lower  attenuation  second 
peak  is  due  to  the  settling  dust  passing  through  the  beam  again.  As  expected,  dust 
attenuation  increases  significantly  with  frequency. 

In  Figure  11  the  angle  dependence  of  the  attenuation  is  examined  for  the  same 
conditions  at  a  frequency  of  45  GHz.  As  expected  the  first  peak  occurs  earliest 
and  the  attenuation  is  highest  for  the  lowest  elevation  angle.  Also  the  separation 
in  time  of  the  peaks  decreases  with  increasing  elevation  angle. 

Jn  Figure  12  the  distance  dependence  of  the  attenuation  is  plotted  for  the  same 
conditions  for  a  fixed  elevation  angle  of  30°.  It  is  seen  that  the  highest  attenuations 
occur  earliest  for  distances  closest  to  the  burst  and  then  decrease  and  occur  at 
later  times  as  the  distance  from  the  burst  is  increased. 
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TIME  (SEC) 


Figure  12.  Attenuation  From  1  MT  Burst  as 
a  Function  of  Distance  From  Burst 

In  general  it  is  seen  that  high  attenuations  occur  at  all  frequencies  as  the  rising 
dust  cloud  passes  through  the  antenna  beam.  For  low  elevation  angles  and  dis¬ 
tances  close  to  the  burst  these  attenuations  may  exceed  100  dB.  Attenuations 
produced  by  the  falling  dust  are  significantly  lower  and  do  not  generally  exceed  5  dB. 

5.  CONCLUSIONS 

Upon  examining  the  effects  of  dust  on  slant  path  propagation  at  millimeter  wave¬ 
lengths  it  has  been  shown  that  the  attenuation  is  heavily  dependent  on  the  maximum 
particle  radius,  the  number  of  large  particles  in  the  distribution  and  the  real  and 
imaginary  components  of  the  index  of  refraction  in  the  range  from  dry  sand  to  clay. 
The  attenuation  is  also  directly  proportional  to  the  fraction  of  the  atmosphere  filled 
with  dust.  Dust  attenuations,  while  low  at  frequencies  below  10  GHz  increase 
significantly  at  millimeter  wavelengths. 

For  the  WESCOM  code,  the  computed  attenuations  behave  as  would  be  expected. 
Very  large  attenuations  occur  immediately  after  the  burst  due  to  fireball  ionization 
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and  dust  when  the  antenna  beam  is  close  to  the  fireball.  At  later  times  the  attenua¬ 
tion  is  caused  mostly  by  dust.  The  attenuation  due  to  oxygen  and  water  vapor  is 
typically  less  than  1  dB  in  the  window  region. 

On  the  basis  of  these  results  it  appears  that  if  the  antenna  beam  intersects  the 
fireball  within  the  first  20  sec  after  the  burst,  then  attenuations  greater  than 
100  dB  are  likely,  and  the  lower  frequencies  are  attenuated  more  than  the  higher 
frequencies  because  the  losses  are  due  predominantly  to  fireball  ionization.  For 
this  reason  it  would  seem  that  the  very  early  dust  model  is  not  too  critical  since 
the  attenuations  are  always  prohibitive.  At  times  from  about  20  sec  to  several 
minutes,  the  period  after  ionization  losses  have  essen+ially  disappeared,  but  during 
which  dust  losses  can  be  tens  of  dB's,  the  accuracy  of  the  dust  model  parameters 
discussed  above  can  significantly  influence  the  results.  Finally,  when  the  dust 
settles,  although  attenuations  of  the  order  of  several  dB  are  likely,  it  would  seem 
that  should  be  possible  to  model  these  more  accurately  than  those  at  earlier  times 
since  the  terminal  velocities  of  the  dust  particles  are  known.  Since  only  the  very 
small  particles  remain  suspended  the  size  distribution  is  less  critical  and  only  the 
density  and  index  of  refraction  of  the  dust  are  important. 
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